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Ql|. Abstract 

■ For a class of random band matrices of band width W, we prove 

regularity of the average spectral measure at scales e > W~ ' 99 , and 
find its asymptotics at these scales. 

in 

1 Introduction 

m 

■<sj- . Define a random operator H on £ 2 (Z) via 

^ . 

fivW^T' u<v<u + W 
H(u,v) = <H(y,u), u>v>u-W (1) 
[0, otherwise , 

^ ■ where the random signs are independent, and the width W G N is a large 



parameter. The integrated density of states iV : R — > [0, 1] is defined by 

N(E ) = lim -( [ ° %(H -E -ie)- 1 (0,0) dE ) , E G R , (2) 

e->0 7T \ y_ 00 / 

where (■) denotes average over the randomness. It is known that the limit 
exists for almost every Eq G R; it is equal to the distribution function of the 
average spectral measure of H (corresponding to the vector Sq.) 
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If iV is (Radon-)differentiable, its derivative is called the density of states, 
and is denoted by p. In this case, 

p(E ) = lim-hiH-Eo-ie)- 1 ^^)); (3) 

the existence of the density of states is equivalent to the existence of the limit, 
and to the boundedness of the expression under the limit. It is believed that 
the density of states exists for any E, is bounded uniformly in W, and admits 
an asymptotic series 

P\Eq) ~ a Q (E Q ) H — I — h • • • , (4) 
where for example 



a (E ) = -Jl-E* . (5) 

7T v 

In fact, there is a natural perturbation expansion using the so-called self- 
energy renormalisation that yields the terms of fll]) one by one; it is equiv- 
alent to the one we describe in Section [2] (see Spencer [11] for the precise 
definitions.) 

However, even the first terms of the series (Tj0) have not yet been rigor- 
ously justified, and it is still unknown whether the density of states exists 
and is bounded uniformly in W — > oo. If the entries of H were replaced with 
random variables with absolutely continuous distribution with bounded den- 
sity, Wegner's estimate would show that the density of states exists; however, 
even in this case, Wegner's argument only yields the width-dependent bound 
p(E ) < CVW. 

We are interested in a simpler problem, namely, the behaviour of the 
expression ({H — E — ie) _1 (0, 0) ) for small e > (depending on W .) 



Theorem. For -1 < E < 1 and e > W~°" , 

a {E)dE 



^(#-£ -z e )-i(0,0))-3 J 



E-Ea-ie 



with a Q as in |3]) ; and C(E ) > independent of W and bounded on any 
interval (— 1 + 5, 1 — 5). In particular, 

^((H-Eo-iey^O)) <C{E ). (7) 
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When e > is fixed and W — > oo, the asymptotics follows from 
the results of Bogachev, Molchanov, and Pastur pj. For e > W -1 ^ 3 , the 
asymptotics (Q follows from the recent result of Erdos and Knowles jl] (see 
also [5] for an extension to more general random band matrices.) These 
results are based on rigorous control of the perturbation series in W. As 
observed by Erdos and Knowles, the perturbative series they use diverges for 
e < W~ 2 ^. 

On the other hand, the methods of [2] and [U [5 J allow to handle the (more 
difficult and physically more interesting) quantity 

(\(H-E -ie)-\0,Q)\ 2 ) , (8) 

and [U [2] actually control the quantum dynamics exp(itH) for t < W 71 / 3 . 

In [7J, Erdos, Yau, and Yin used different methods to prove an analogue 
of (|6]) for finite band matrices of size N x N, for e > W~ l \og c N, and with 
error term < (We)~ 1 ^ 2 \og c N. They also control the fluctuations of 

(H-E Q -ie)- l {Q,Q) 

at these scales. 

In this note we go back to the perturbation series (jl]), and introduce a 
regularisation procedure which allows to justify it at scales e > W^ -0 - 99 . For 
now, we only deal with the average spectral measure; however, we believe 
that the method could also be applicable to OH]). One advantage of the 
perturbative method is that it allows to get the (optimal) error term 0(1/1^); 
actually, our method can be used to justify the first W om terms in the 
asymptotic expansion, and thus obtain an approximation with error term 
which is exponentially small in the width W . 

We refer to the survey of Spencer [TT] for other results and problems 
pertaining to random band matrices. 

The paper is organised as follows. In Section |5] we describe the formal per- 
turbation series. In Section [3] we collect several technical statements which 
are used in the proof. Section |4] briefly summarises the topological classifi- 
cation of paths, based on pj. In Section |5] we estimate the contribution of 
every equivalence class to the perturbation series. The proof of the theorem 
appears in Section El and is followed by the concluding remarks of Section [7] 

Acknowledgment. Tom Spencer encouraged me to study the average spec- 
tral measure on short scales, suggested several crucial steps in the argument, 
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and commented on a preliminary version of this paper. I thank him very 
much. 



2 Formal perturbation series, and combina- 
torial preliminaries 

Denote by 1_{W) the graph having Z as its set of vertices, and u ~ v if 
< \u-v\ < W. 

Let T n and U n denote the Chebyshev polynomials of the first and second 
kind, respectively: 

t t a\ ( a\ tt < a\ sin (( ra + X ) g ) / Q \ 
l n (coso) = cos(nU) , U n {cosv) = : — ; (9) 

sin u 

we formally set U- 2 = U-\ = 0. Let 

U n , w = U n - ^l—^J n _ % . (10) 

Then (see [9l Lemma 2.7]) 

U n , w (H)(u ,u n ) = (2W -l)- n/2 J2 H ( u o,ui)H( Ul ,u 2 )---H(u n ^,u n ) , 

where the sum is over (n + l)-tuples (u , U\, u 2 , ■ ■ ■ , u n ) such that Uj ^ Uj +2 , 
j = 0,1,2, ■■■ ,n — 2 (non-backtracking paths.) In particular, 

(C/ niW ,(F)( Wo ,«n)) = (2W - l)-"/ 2 Paths„K, M „) , (11) 

where Paths n (wo, u n ) is the number of paths from uq to u n in Z(W) which 
pass every edge an even number of times and satisfy the condition above (in 
particular, Paths n (u, v) = 5 UV .) Combining ffT0|) with ffTTj) . we obtain: 

n/2 

(U n (H)(u Q ,u n )) = (2W - ir n/2 ^ Paths n _ 2m ( Mo , Mn ) , (12) 

m=0 



and hence 



n/2 

(U n (H)(u ,u )} = (2^-l)- n / 2 ^Paths n _ 2m K,«o) • (13) 



m=0 



Observe that Paths m (-u , u o) does not depend on u G Z, therefore we denote 
it simply Paths m . Also, Paths m = if m is odd. 
Next, T n = (U n — U n -2)/2, therefore 

(T n (tf)(0,0)> 

«/2 . (14) 



2{2W 



\ ^ n/2 ^ |Paths„_ 2m - (2W - 1) Paths n _ 2m - 2 | , 



m=0 

where we formally set Paths m = for m < 0. 

It will be convenient to rewrite |Paths„_ 2m — (2W — 1) Paths„_ 2m _ 2 | 

in a different form. Let Paths" be the number of paths which pass every 
edge an even number of times and satisfy the strengthened non-backtracking 
condition u ^ w 2 , u\ ^ u 3 , . . . , u n -\ ^ U\. Then 

Paths n = Paths" + (2W - 2) ( 2W ~ 1 ) i_lpaths °-2 J - • 

l<j<n/2 

Therefore 

Paths n - (2W - l)Paths„_ 2 = Paths" - Paths°_ 2 . (15) 

To make this identity valid for all n > 1, we formally set Paths^ = 2W — 1, 
Paths" 1 = 0. 

Now, 5(E — Eq) admits a formal expansion 

5{E - E ) 7 =L= { 1 + 2 Tn(E )T n (E) 1 (16) 

We shall discuss a regularised version of this series in the next section; for 
now, we remark that (f!6l) is a rigorous identity in L 2 (— 1, 1), since T n are the 
orthogonal polynomials with respect to the measure — f E . 

Convolving ( f!6l) with p, we obtain: 
p(E ) ~ (5(H - E )(0,0)) 

^( 1+ #' (£JWi ' ,W » ) J (17) 



oo 



_± = ) 1 + y T 2n (E Q) y r Paths o _ Paths o i 



7T 
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Picking the addends which are not divided by powers of 2W — 1, we see that 
the leading term is 



nyfl^Ei I 2W-1 v J 7i 



x (-2W + 1) > = -\ll-El 



In this way one can also obtain the full expansion of the form (J3J) which, 
unfortunately, diverges. 

One may introduce a regularisation factor and consider the expression 

-{%{H-E o -ie)- l (0,0)) 

71 

7 =i=(l + 2f^e--T n (E )(T„(iJ)(0,0))l , (19) 

^ v o I ^ J 

which corresponds to the density of states averaged over an interval of width 
~ e about Eq (more precisely, the convolution of the density of states with an 
approximate ^-function of width e.) However, it also diverges for sufficiently 
small e. The reason is the large contribution to (T n (H)) of the part of the 
spectrum of H outside [—1, 1]. Namely, it is known [10J that a part of the 
spectrum lies at distance ~ W~ A ^ from [—1, 1]; the polynomials T n grow as 



\T n (±(l + 5))\ ^exp(V25n) , 

therefore the series above can not converge for e <C W~ 2 ^. 
Here we propose a different regularisation: 



1 + 2 J2 v(ne)T n (E ) (T n {H)(0, 0)) 



n=l 



and justify its convergence for e > IV -0 - 99 . For reasonable <ft, we shall see 
(in Lemma I3.ip that it also represents an average of the density of states 
over an interval of width f=s e about Eq. Therefore a posteriori we obtain an 
expansion for 

-{^(H-Eo-ier^O)} . 

The main step is to pick (p carefully. In particular, the argument above 
shows that (p has to decay faster than exponentially at infinity. As we shall 
see, it will also be convenient to have (p positive-definite and analytic. 



6 



3 Some auxiliary statements 

Let ip : IR — > R + be a smooth even function such that ip(0) = 1, ip decreases 
to zero on R + , and tp > 0. We shall take 



1 f°° 

p q (t) = — ew [-s^-(t-s) 2 "}ds, 

A q J-oo 



(20) 



where q is a (large) integer parameter which we shall choose later, but the 
next lemma will be applicable in the general setting. 
For e > and Eq = cos 9 £ (—1) 1)) set 

f E(h€ (E) = 1 + 2J2 tp(ne)T n (E ) T n (E) . 

n>l 

Dirichlet's criterion shows that the series converges. 
Lemma 3.1. 



/cos o ,e (COS 0) 

1 oo 



2e ^ |"U 

m=— oo 



m 



+ Oq 
2tt 



' > 7 



m 



0-Oq 
2tt 



• (21) 



Note that, if <p > 0, the right-hand side is an approximate 5-function in 
9 of width ps e. 



Proof. First, 



/COS 00,6 (£) = ~ ^ #(™) , 



where 



Now, 



$(t) = y?(te){ exp(it(# + 0„)) + exp(it(# - O ))} • 



/oo 
<?(*) exp(-2mt£)dt 
-oo 



i r . n 



2tt 



0-0q 
2vr 



(22) 



therefore the lemma follows from Poisson's summation formula. 



□ 
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Now we need some properties that are specific for (p q from (I2(jp . Denote 

/oo 
exp(-2nix£ - x 2q )dx . (23) 
-oo 

This is obviously an entire function. The following lemma can be proved 
using a saddle-point argument: 

Lemma 3.2. For any 5 > there exist Cs,c$ > such that 

\F q (pe l *)\<Csexp{-csp^} 
for ~f q + 5 < <P < f q - 5 and for n - f q + 5 < < vr + f q - 5. 



Now take <p = <p g in Lemma I3.lt then ip q = F q . From Lemma I3.2[ 
both sides of ( 12 ip are analytic functions of 8, therefore the equality can be 
extended to complex 9. In particular, f Eo>e is bounded (uniformly in e > 0) 
on an open interval containing [—1, 1]. 

Finally, we state - for use in Section [5] - the definition and a couple of 
properties of divided differences; we refer to the survey of de Boor [3] for the 
proofs. 

For (distinct) Zi, ■ ■ ■ ,ze in the domain of definition of a function /, con- 
sider the (E — l)-th divided difference f[z\, • • • , ze], defined as follows: 

f[zi]=f{zi), f[zi,---,z E ] = • 

z\ — ze 

An equivalent definition is given by 



f[zi,--- ,z E } = Y^ 



Lemma 3.3. Let f be an (E — 1) -times differentiable function in a convex 
domain D C C. Then for any distinct Zi,-- - ,Ze € D there exists z* G 
conv(zi, • ■ ■ , ze) such that 

,r , / (E - 1) (^) 

J[zu--- i z E\ 



(E-l)\ ■ 



8 



Lemma 3.4. For z\, ■ ■ ■ , Zg E C, 

E E 

(ni,-,n E )6A Bi „ e=l e=l 

w /i ere 

and m n -i(z) = z n ~ x . 



| Til, • • • ,TIe > 1 | Tli + h Tl E = Tlj 



4 Classification of paths 

Consider the collection Paths^ of paths u = 0,ui,u 2 , ■ ■ ■ ,u 2n -i,U2n = 
in Z(W) such that every edge appears an even number of times, Uj ^ Uj +2 
for j — 0, 1, ■ • • , n — 2, and w n _i ^ u\. 

A pairing of a path is a paring of 0, 1, • - • , 2n — 1, so that if j is paired 
to f , then either u^- = Uj> and = iZj/ +1 , or Uj = Uj> + i and Uj +1 = Uj>. 
Every path in Paths^ has at least one pairing. 

Now we divide the couples (path, pairing) into equivalence classes. The 
procedure is a slight elaboration of the one from [HI Section II. 1]; it works as 
follows. 

We look for a pair (j, j') so that j is paired to j 1 and j + 1 is paired to 
either j 1 + 1 or j' — 1, and unite the j-th and the j + 1-th edge into a single 
one (and the same for their counterparts.) Continuing this process, we arrive 
at a multigraph G = (V, E) with a marked vertex Vq (corresponding to 0) 
together with a path p which passes every one of its edges twice. We call the 
equivalence class D — (G — (V, E),p) a diagram of order 1. 

Paths which contain edges passed more than twice admit more than one 
pairing, and therefore correspond to more than one diagram. If a path passes 
a certain edge 4 times, we correspond to it a diagram of order 2 in a similar 
way, and so forth. Thus, every path which corresponds to a high-order 
diagram also corresponds to a sequence of diagrams of lower order, and the 
total number of paths in Paths2 n can be computed using the inclusion- 
exclusion formula: 

Paths 2n 

= ^^( — iy +1 # {paths of length 2n corresponding to D} . (24) 

j>l D of order j 
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Let us compute the number of paths corresponding to a diagram D = (G,p). 
Fix an arbitrary ordering of the edges in E. To construct a path correspond- 
ing to D in Z(W), we first choose the positions R v G Z corresponding to 
v e V, and the number n e of edges corresponding to every e G -E, so that 

R vo = and % + h He — n. 

Then, for e = (u,v) from 1 to E, we start a random walk from R u 

conditioned not to backtrack, and denote by P e _R„ j the probability 

that it arrives at R v after n e steps without violating the non-backtracking 
conditions created by the edges corresponding to / < e from the previous 
steps at its first and last step. The number of paths corresponding to D is 
then equal to 

E 

2W(2W-l) n - 1 ^2^2l[Pe{Ru^Ru} ■ 

R n e=l 

One can proceed using the expression due to Smilansky [T2] for P e in 
terms of the transition matrix P of the random walk on ~L(W). Instead, let 

us denote by P ^R u ~% -R„ j the probability that a random walk from R u and 

conditioned not to backtrack arrives at R v after n e steps. Then 

P e {Ru^Rv}^P{Ru^Rv) ■ 

More formally, the left-hand side of the last equality can be expressed as the 
right-hand side plus a sum of similar terms of the same form with different 
parameters. Therefore we can essentially regard the above approximation as 
an identity. 

The degree of the marked vertex in a diagram of order one is at least two, 
and the degree of every other vertex is at least three. If these inequalities 
are saturated, the diagram is called simple. For simplicity, let us focus on 
simple diagrams of order one (which were called diagrams in [8].) 

The genus of a diagram is defined as j(D) — E — V + 1. A diagram of 
genus 7 satisfies E = 37 — 2, V — 27 — 1; the number .0(7) of diagrams of 
genus 7 satisfies (see [S]) 

(7/C) 7 < D(i) < (C7) 7 . 

We remark that there is exactly one diagram of genus 7 = 1; it contains one 
vertex, and one edge (which is a loop connecting the vertex to itself.) 
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5 Embeddings into Z(W) 



Fix a multigraph G = (V, E) with a marked vertex Vq G V. For g G C with 

|gr| = 1, set 



Emb(G) = Emb(G;#,e) 

= E E ^(E^ e ) II ^ e Pe{^-^} 

R n:_B^N eS-B e=(u,v)^E 

where the exterior sum is over R : V — > Z such that R(vo) = 0. 
Our goal in this section is to prove 

Proposition 5.1. 

1. For any g with \g\ = 1, 



(25) 



where C(q) > is a constant depending only on q. 

2. If G is the multigraph corresponding to the (unique) diagram of genus 
7=1, the same bound holds without the logarithmic factor. 

Remark. The logarithmic factor is probably unnecessary in the general case 
as well; this is however not essential for our purposes. 

To make the computations more transparent, we make several simplifica- 
tions. First, set 

Emb(G) = Emb(G;#,e) 

R n:_B->N eS-E e=(u,v)&E 

By the argument sketched in the previous section, it is sufficient to prove 
the bound for Emb. 

Let P be the transition matrix of the usual random walk on 1_{W). From 

n. 

p {*» 54 M = v(2w-ir "- w ( vfe) {R - • (27) 
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where the last brackets stand for taking matrix elements. Instead of Emb(G), 
we shall prove the bound for 

Emb # (G) = Emb*(G;g,e) 

=E E ^(E^) n g ne P ne (Ru,R v ) ■ ^ 

R n:£->N egE e=(u,v)£E 

Using (j2"T|) . one can repeat the argument and obtain the same bound for 
Emb(G) (and hence also for Emb(G).) 

We start with a representation of Emb*(G) in Fourier space. The oper- 
ator P is translation- invariant, therefore diagonal in Fourier space: setting 
e^(n) = exp(27ri£n), < £ < 1, we have: 

Pe t = w({)e £ , (29) 

where 

w ® = W E cos(27tjO = cos(vr(^ + 1)0 . 

W ^ W sin (7rt ) 

3=1 

For future reference, we remark that 

£ l + c Iym in K ,l-0 ' < 3 °) 

Choose an ordering of the vertices of G, and for every edge e = (u, v) G E, 
u -< v, introduce a variable £ e = if u ^ v, set £( U)U ) = -£(«,„). 

Lemma 5.2. S'et S^z) = 2~^ n >i <f q (ne)z n ~ 1 . JTien 

m # (g) = /•••/ M^w^r-^y n^fe)) > 



[0,1] 



where 6 Kirch is the Lebesgue measure restricted to the (E — V + 1)- dimensional 
subspace defined by the Kirchhoff constraints 

Vw g v £(«.«) = • 
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Proof. From 

P= I w{Oe^e^ , 
Jo 

and ^ 

P n (R 1 ,R 2 )= [ w(O n exp(2vr^( J Ri- J R 2 ))^ • (31) 
Jo 

Now substitute ([SI]) into (J2SJ. We obtain: 
Emb # (G') 

= EE^E ne£ ) /■■•/ II [^ e (^fe)) ne exp(27rze e (i? u -i?,))" • 

R n e j Q ^p e=(w,«) 

Exchanging the summation over R with the integral, we see that 

Emb#(G) =X>(5>"0 /■■■/ ^Kirch(e) II (^fe))™ 6 • ( 32 ) 

n e [0,1] B e=(«,«) 

Now we can exchange the sum over n with the integral. According to 
Lemma 13. 4[ 

e n^^r = • • • ,gw(t E )] n^fe)) > 

niH \-riE=n e e 

therefore 

n e e=(u,v) 

n>l e 



□ 



According to the mean- value theorem (Lemma 13. 3p . 



Se[gw&), ■ ■ ■ ,gw(£ E )} = (33) 

for some — 1 < w* < 1. To conclude the proof of Proposition 15. 1[ we need 
one more lemma: 
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Lemma 5.3. For j > and \z\ < 1 , 

where the constant C(q) does not depend on e > 0. 

Proof of Proposition 15. 1\ Applying Lemma [5.21 the relation (J33l) . and Lem- 
ma HQ} we obtain: 

|Emb # (G)| < {^J^j E+1 J «)H HCeM ■ (34) 

The integration is over an (E — V + l)-dimensional subspace, therefore the 
bound fl30|) concludes the proof of 1. 

The estimate 2. can be verified directly from Lemma 15.21 

□ 



Proof of Lemma \5. 3[ First, 



/oo /*oo 
exp(27m;£) <p q (£) d£ = / exp(2mneg) F q (g) 
■oo J — oo 



dt 



Fix < <p < tt/2 satisfying the assumptions of Lemma [3.21 and deform the 
contour of integration to L = | arg£ E {(j), it — 0}}. Now 

S e (z)= [ F g (0 2 ^exp(2 7 rm60^- 1 ^ 

Jl n>l Aq 

f „ (n 2 exp(27rieQ 
Jl ^ (l-zexp(2me0) A q ' 

and 

2 exp(27rie£) d£ 



S?(z)=j\ / F q (0 



L (1 - ^exp(27rkO) i+1 A 



Taking absolute values and applying Lemma I3T2| we conclude the proof. 

□ 
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6 Proof of Theorem 

Let q be a large natural number, and let 77 > be a small real number; we 
shall choose them later. We shall work with the function ip q from (1201 . or 
rather with its truncated version 

<p(t) = <p q (t)l{\t\<wn}. (35) 

Denote 

f{E) = 1 + 2J2 V q (ne)T n (E Q ) T n (E) . (36) 

n>l 

Lemma 6.1. For any q > 50 and 77 < 1/100 

(7(10(0. 0)) = 1 + <p q (2e)(l - 2El) + 0{l/W) . (37) 

Proof. First, 

w — 1 n ° 

= 1 - Vq (2e)(2E 2 - 1) 2 + 2^^(2ne)T 2n (^)(T 2B (^)(0,0)> , 

' n=2 

where n = L^V e J- B Y dH, 

(T 2n (g)(0,0)) = ^^yp Paths° B + [•••" , 

where the last brackets enclose the terms of higher order which can be anal- 
ysed similarly to the leading term, Now apply the classification of paths 
described in Section [2j To simplify the notation, we explicitly write the con- 
tribution of simple diagrams, and collect all the rest in the remainder term. 
This yields: 

Paths°„/(2W<--1)" 

= E E E e n 

7<no 7(D)=7 R nxH Vn E =n e=(u,v) 



(38) 
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therefore, 



(2W- l) n 



7<no 7(D)=7 n<rao 

E E II h{R,-^} + 

R niH hn£=n e =(u,») 

It is not hard to see (cf. [10]) that 



(39) 



R mH hn B =n e=(u,«) 



< Cexp [-cn 2q e 2q ] 



(Cn 



57-4 



' 57-4 1 



(40) 



therefore the left-hand side of ff4U|) is very small for n > n if 

2g- 1 



2? 



> 0.99 . 



(41) 



Under this condition the sum over n in ( 13 9 p can be extended to infinity. Thus 

(7(i0(o,o)> 

= l-<^(2e)(2£ 2 -l) + Emb(G) +[•••] +0(1/W) . 

7<n D=(G,p),7(D)=7 

For D of order 7, Proposition 15.11 and the subsequent remarks yield 



Emb(G) 



< 



C(q) 



3-1 



logVT 



(42) 



and the logarithmic factor is redundant for 7 = 1. As n = [W v /e\ <C W 
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the estimate (1421) implies: 

(/(tf)(0,0))-l + ^(2e)(2£ 2 -l) 



< 



C 



\l-g\W 
0(l/W) . 



+ £ (^) 7 (irr5i) 1, iF + ^» (43) 



2<7<ng 



□ 



The Chebyshev polynomials T n satisfy 



max \TJE)\ < (CW) n/2 

Be [- W l/2 iff l/2] 



Therefore for -W 1 ' 2 < E < W 1 ' 2 



<p g {ne)T n {E )T n {E) 



n>no 



< Cexp [-cn 2q e 2q + Cn\nW] = 0(1 /W) 

n>no 



as long as 2qi] > r\ + 0.99, or: 



The spectrum of H lies in 

W 



q > 



r/ + 0.99 



2r] 



(44) 



W 



c [-iy 1/2 , 



/2l 



1 - 1 ' V2W-1 

thus, under the assumption f T4"4"|) . the conclusion of Lemma T6. II remains valid 
for f Eo>e in place of /. 

Next, one can replace fE 0j6 with /s ,e|r 1 ^ (using Lemma I3T2"} the remark 



following it, and the fact [TO] that the density of states is small outside 
[— 1, 1].) Since 



E o,qr_i 



T-£ 2 )>0 



is an approximate 5-function of width e at E , and E , e are arbitrary (subject 
to the constraint e > W~ ' 99 ), we can replace it with any other approximate 
5-function h Eoe , such as the Stieltjes kernel that appears in the statement 
of the theorem. Indeed, for any e = ^(W) ^> e(W), and any — 1 < E < 1, 
one can approximate J h^ Q ~dN by positive linear combinations of J h Eo ,t dN 
(with different Eq.) □ 
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7 Two remarks 



i. As we remarked in Section [2J the divergence of the "naive" perturbation 
series (obtained from the self-energy renormalisation procedure) follows from 
the divergent contribution of the spectral edges. It is probable that a similar 
reason is responsible for the divergence of perturbation series also in other 
problems, such as the density of states in the Anderson model (see Erdos- 
Salmhofer-Yau [6].) 

ii. The restriction e > W~ ' 99 in the main theorem appears for the following 
reason. It is an artefact of the approach that ()6]) can be justified for a given 
e > only together with the first wl/e terms of (j3J). However, only the first 
~ W terms of fll]) are reasonably small (say, smaller than 1), therefore we do 
not see how to make the current approach work for e <C W~ x . It is possible 
that the power —0.99 can be improved to —1 using a more careful choice of 
the test function ip. 
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